Transition between Quantum States in a Parallel- Coupled Double-Quantum-Dot 



J.C. Chen 1 , A.M. Chang 1 and M.R. Mclloch 2 

1 Department of Physics, Purdue University, 
West Lafayette, IN 47907 
2 School of Electrical and Computer Engineering, 
Purdue University, West Lafayette, IN 47907 

(Dated: February 6, 2008) 

Strong electron and spin correlations in a double-quantum-dot (DQD) can give rise to differ- 
ent quantum states. We observe a continuous transition from a Kondo state exhibiting a single- 
peak Kondo resonance to another exhibiting a double peak by increasing the inter-dot-coupling (t) 
in a parallel-coupled DQD. The transition into the double-peak state provides evidence for spin- 
entanglement between the excess-electron on each dot. Toward the transition, the peak splitting 
merges and becomes substantially smaller than t because of strong Coulomb effects. Our device 
tunability bodes well for future quantum computation applications. 
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The double-quantum-dot (DQD) is emerging as a ver- 
satile system for studying a variety of strongly correlated 
behaviors OIQEIEIIHSB]- Following the experimental 
demonstration of the Kondo imp urity-sp in screening ef- 
fect in single quantum dots pi. lij. ITol lTll H2I IT^|. recent 
theoretical investigations of the coupled-DQD system is 
uncovering new correlated behaviors 0, y, 0, IE • 
These works suggest that the DQD enables a realization 
of the two-impurity Kondo problem first discussed in the 
context of metallic systems 0, [Til LLj- LLJj m which a com- 
petition between Kondo correlations and antifcrromag- 
netic (AF) impurity-spin correlation leads to a quantum 
critical phenomenon. In a different regime of parameters, 
a related quantum critical phenomenon can occur driven 
by a competition between intra-dot Kondo coupling to 
leads and the inter-dot-coupling 0,0. In each scenario, a 
transition is predicted to occur between a quantum state 
characterized by a single-peaked Kondo resonance, and 
a different state with a double-peaked resonance. De- 
pending on model and DQD gcomctry-whcthcr series or 
parallel coupled-both a continuous or discontinuous 
behavior in the Kondo peak characteristics have been 
predicted. The quantum transition in the two-impurity 
Kondo problem has received wide attention in the theo- 
retical literature in the past two decades, to a large ex- 
tent because the Kondo to antifcrromagnctic transition 
involves an unusual non-Fermi liquid fixed point. Ex- 
perimental investigation of this problem thus far has not 
been reported. 

Here we describe transport properties of an artificial 
molecule formed by two-path, parallel-coupled doublc- 
quantum-dots, where the inter-dot-tunncl-coupling, t, 
can be tuned. In the Kondo regime the differential- 
conductance, dl/dV, exhibits a single peak centered at 
zero-bias for t comparable to the lead-coupling induced 
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level broadening. Increasing t by less than 10% resulted 
in a continuous evolution into a split Kondo resonance. 
At the same time the conductance at zero-bias exhibits 
a maximum in the vicinity of the transition. This peak 
splitting behavior in conjunction with distinct tempera- 
ture dependences in the different regimes demonstrates a 
direct observation of an inter-dot-coupling-induccd quan- 
tum transition. Moreover, on the double peak side the 
zero-bias conductance becomes suppressed; this suppres- 
sion represents direct evidence that the localized dot spins 
are becoming entangled into a spin singlet. 

While a double-peaked, coherent Kondo effect was 
observed by Jeong et al. in a series-coupled DQD 
geometry [T3, the existence of a singlc-peakcd Kondo ef- 
fect at finite inter-dot-coupling has not previously been 
established. The series-geometry is unsuitable for an in- 
vestigation of the transition because even a slight de- 
crease in t can drastically reduce the tunnel current be- 
low detection. Therefore neither a single-peaked behavior 
nor a transition could be observed. Furthermore, direct 
evidence for spin-entanglement was not obtainable. In 
the parallel geometry the undesirable suppression of the 
conductance is avoided. 

The Kondo effect in a single QD results from the cou- 
pling between the dot excess (unpaired) spin and the 
spin of the conduction electrons in the leads. The en- 
ergy scale is given by the Kondo temperature Tk ~ 
VUTexpl— 7r|/x — e \(U + eo)/r?7] where U is the on-site 
charging energy, eo is the energy of the single-particle 
level, T reflects the dot level broadening from coupling 
to the leads, and /i is the chemical potential. The fully 
symmetric DQD system contains the additional inter- 
dot-coupling parameter, t. The two magnetic impuri- 
ties, realized by a single excess-spin on each dot, interact 
through an effective AF coupling, J = 4t 2 /U. The new 
energy scales, t and J, together with U, eo, and T, in- 
troduce a rich variety of correlated physics. When U is 
the dominant energy scale as in the case of our quan- 
tum dots, for energies below U two different scenarios 
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arc predicted for dots coupled in series Jl], 0, 0, LH El- 

When t< T, the system can be cast [l| into the two- 
impurit y K ondo problem initially discussed by Jones et 
al. [l4L \13L jig. The competition between Kondo effect 
and antifcrromagnetism appears as a continuous phase 
transition (or crossover) at a critical value of the cou- 
pling J/Tk ~ 2.5. If however, t is tuned to t > T before 
the AF transition point can be reached, the system un- 
dergoes a continuous transition from a separate Kondo 
state of individual spins on each dot (atomic-like) to a co- 
herent bonding-antibonding superposition of the many- 
body Kondo states of the dots (molecular-like) 0, Li • 
Both the AF state and coherent bonding state exhibit a 
double-peaked Kondo resonance in the differential con- 
ductance versus source-drain bias and involve entangle- 
ment of the dot spins into a spin-singlet. Therefore they 
are likely closely related to one another. The parallel- 
coupled case has only recently been analyzed in a model 
without inter-dot-tunnel-coupling, t, where the AF cou- 
pling occurs via electrostatic coupling. A discontinuous 
jump in the differential conductance is predicted when 
the AF state becomes favored @- 

Our device was fabricated on GaAs / Al x Ga\- x As het- 
erostructure containing a two- dimensional (2D) electron 
gas 80nm below the surface, with electron density and 



mobility of n = 3.8 x 10 1 



and fj, = 9x I0 5 cm 2 /Vs, 



respectively The lithographic size of each dot is 170nmx 
200nm (see Fig. 1(a)). The eight separate metallic gates 
are configured and operated with five independently tun- 
able gate voltages. The dark regions surrounding (and 
underneath) gates 5 represent 120nm thick over-exposed 
PMMA, which serve as spacer layers to decrease the lo- 
cal capacitance thus preventing depletion, enabling each 
lead to simultaneously connect to both dots. The experi- 
ment carried out at a lattice-temperature of 30mK. Stan- 
dard, separate characterization of each dot in the closed 
dot regime [13 yielded a charging energy U of 2.517meV 
(2.95mcV) for the left (right) dot, with corresponding 
dot-environment capacitance Cs ~ 63.6a/(54.3a/), and 
level spacing AE w 219^eV(308fieV). Modeling the dot 
as a metal disk embedded in a dielectric resulted in disk 
of radius r e « 70nm(60nm) with 58 (43) electrons. To 
characterize the DQD and demonstrate the full tunability 
of our device, in Fig. 1(b) we show the Coulomb blockade 
(CB) charging dia gram of the conductance versus plunger 
gates V2 and V4 [lg as t was increased, for weak cou- 
pling to the leads where Kondo correlation is unimpor- 
tant. The central pincher gate V5 controlled t where a re- 
duction of the gate voltage decreased t. At weak coupling 
(t small) , the electrons separately tunnel through the two 
nearly independent dots forming grid like pattern, yield- 
ing rectangular domains in Fig. 1(b)(1). With increasing 
t charge quantization in individual dots is gradually lost 
as the domain vertices separated and the rectangles de- 
formed into rounded hexagon. At large t when separate 
charge quantization is fully relaxed, the two dots merge 
into one and the domain boundaries become straight lines 
(Fig. 1(b)(6)). The evolution of the conductance pat- 
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FIG. 1: Device and characterization in the CB regime: (a) 
Scanning electron micrograph of the device, (b) Logarithm of 
double dot conductance as a function of gate voltages V2 and 
V4 for weak lead-dot coupling. The color scale indicates the 
magnitude of conductance. The voltages in central pincher 
gate V5 are (1)-0.7477V (2)-0.6573V (3)-0.6534V (4)-0.6504V 
(5)-0.6494V (6)-0.5701V. 



tern demonstrates the tunability of our DQD from ionic 
to covalent bonding states and bodes well for quantum 
computation applications |l9j |. 

The parameters T and t govern the delicate DQD 
Kondo physics. To obtain the necessary large t, the cen- 
ter pincher gate V5 was set so that the zigzag pattern in 
the charging diagram is barely visible (Fig. 2(a)), ensur- 
ing that the Kondo valleys can be located. To accom- 
plish the formation of Kondo states in both dots, pincher 
gates VI and V3 were tuned to give a sizable T w AE 
to ensure strong Kondo correlation. An estimate for t is 
based on the fact that the charging diagram indicates a 
configuration close to the limit of a merged, single large 
dot (Fig. 2(a), so that the level broadening TT\t\ 2 /AE 
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FIG. 2: Device characterization in the Kondo regime: (a) 
Charging diagram for the third cool-down depicted in a gray- 
scale plot of the conductance crest as a function of gate volt- 
ages V2 and V4 (V5= -0.5965V). (b) Differential conductance 
traces from valleys 1 to 6 in Fig. 2(a). The insets indicate 
the spin configuration of the uppermost, occupied electronic 
levels on each dot. Note that a single upward-pointing ar- 
row only denotes an unpaired electron, and is not intended to 
represent the actual direction of spin alignment. 



should be comparable to the level spacing AE, yielding 
t w 150/j,eV, while T is deduced to be ~ l50/j,eV from 
the half-width of the Coulomb blockade peaks. 

Care is required when changing t. The mutual capac- 
itive coupling between the gates and dots gives rise to 
a complex capacitance matrix. Operating V5 to tune t 
simultaneously affected the charge on the dots and other 
gates. Furthermore, slight residual asymmetry in the 
coupling to the leads caused the maximum of the Kondo 
zero-bias anomaly (ZBA) in the dl/dV to occasionally oc- 
cur at nonzero voltages Vsd (termed the anomalous 
Kondo effect |2(j ). By adjusting a combination of the 
remaining gates we tuned the ZBA to be nearly symmet- 
rical about zero-bias. For practical purposes we relied on 
plunger-gates V2 or V4, which were found to follow V5 in 
an approximately linear manner. This procedure can be 
thought as experimentally diagonalizing the capacitance 
matrix. The desired configuration of an unpaired excess- 
spin on each dot could be maintained within a tuning 
range of ~ 4 ~ 7mV in V5 without causing a sudden 
change in the charge configuration. Note that V5 is set 
typically sw 70mV above pinch-off. Therefore the relative 
tuning range in this Kondo regime is roughly 6-10% of 
closure and even smaller for the corresponding fractional 
change in t, since tunnel-coupling becomes exponentially 
suppressed in the small value limit. 

Three distinct spin configurations may appear assum- 
ing even-odd electron filling and focusing on the top- 
most states in each dot (Fig. 2(b)). Our investigation 
was carried out in such regimes. Non even-odd behavior 
was also observable, but will not be discussed |2l|. The 
Kondo valleys and spin states in each dot were identified 
by measuring the differential conductance, dl/dV, versus 
source-drain bias, Vsd, at an electronic temperature of 
~ AOmK. For regions 1-6 in Fig. 2(a), we observed the 
expected appearance and disappearance of Kondo reso- 
nance peaks near zero-bias shown in Fig. 2(b). 
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FIG. 3: The differential conductance, dl/dV, versus Vsd 
for different inter-dot coupling strength tuned via V5 gate, 
(a) From top to bottom: V5= -0.6155, -0.616, -0.6165, - 
0.617, -0.6175, -0.618, -0.6185, -0.619, -0.6195, -0.62V, second 
cool-down, (b) From top to bottom: V5= -0.5965, -0.5970, 
-0.5975, -0.5978, -0.5982, -0.5987, -0.5990, -0.5992, -0.5995V, 
first cool-down, offset by 0.005e2/h for each trace, (c) From 
top to bottom: V5=-0.5925, -0.5928, -0.5931, -0.5934, -0.5937, 
-0.5940, -0.5943, -0.5946, -0.5949, -0.5952, -0.5955, -0.5958, 
-0.5961, -0.5964, -0.5967, -0.5970, -0.5973, -0.5976, -0.5980, - 
0.5983, -0.5986V, third cool-down as measured in Kondo val- 
ley 2 of Fig. 2(b), offset by 0.02e 2 /h for each trace. The 
traces vary in different cool-downs due to the rearrangement 
in the occupation of charged in defects and traps, (d) Selected 
data from (c) without offset. Bottom half shows first quan- 
tum state regime (from top to bottom: V5=-0.5961, -0.5964, 
-0.5970, -0.5986). The upper half shows the second regime 
(from top to bottom: V5=-0.5946, -0.5940, -0.5937, -0.5934). 
(e) ZBA peak amplitude g = dI/dV\v=o, peak splitting SV 
and width W extracted by first subtracting three types of 
background signals followed by a two- Gaussian fit. The fit 
and data are virtually indistinguishable. The simulated back- 
grounds are •; two exponential functions + constant (e.g. 
marked by the two gray dash lines for the bottom trace in 
(c) and combined functions shown as gray dot line); A two 
Boltzmann and ■ linear function. The bottom figure shows 
the two peak widths obtained after the subtraction of a two 
exponential background. Results for a single Gaussian fit (V5 
< -0.594V only) are indicated by open-circles (o). Typically, 
the fitting error is smaller than the symbol. 
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When t was tuned two distinct regimes of behavior 
were evident in dl/dV (see Fig. 3(a)-(d)). The main fea- 
tures in the first regime (see Fig. 3(a), (d) bottom half) 
were the clear presence of a single peak in the ZBA, an 
increasing peak width and an increasing linear conduc- 
tance g = dl / dV\v=o with increasing t. In the second 
regime where t was increased further (see Fig. 3(b), (d) 
upper half), the single ZBA peak developed into two 
peaks. In some Kondo valleys, both types of behavior 
were observed as shown in Fig. 3(c),(d), in which the 
transition is seen to take place in a continuous manner. 
In the transition region, the ZBA peak broadened and 
its zero-bias value, g, approached a maximum, became 
flat before dropping as the ZBA peak gradually split into 
two. The suppression of g on the double-peaked side can 
be attributed to the emergence of spin-singlet correlation 
between the two dot spins 0, H H (l E UM Thc maxi- 
mum g is about ~ 0.1e 2 /h, reduced below the theoretical 
unitary limit of Ae 2 /h. This reduction can occur when 
the energy levels in the two dots differ and the dot-lead 
coupling is asymmetric 0, H3 • 

Analyzing the transition with detailed curve fitting 
we quantify the amplitude, splitting and width of the 
ZBA peak(s) at each V5 gate voltage within well-defined 
bounds. In the absence of precise theoretical functional 
forms either for the conductance background or the ZBA 
peak itself, we first systematically subtracted three sen- 
sible backgrounds. After subtraction of each background 
we fitted the ZBA using three (combined) peak shapes: 
one (or two) Gaussian, Lorentzian and the Breit-Wigncr 
resonance form. Thc resultant parameters invariably ex- 
hibited similar trends as shown in Fig. 3(e). Thc double- 
peak feature visibly disappeared at V5 ~ —0.5941^ as 
indicated by the dash line, and roughly coincides with 
the g maximum position. Note that the peak splitting, 5, 
dramatically reduces from a maximum of ~ 120fieV to 
~ when V5 is slightly reduced from -0.5925 V to ~ 
-0.595 and correspondingly t is reduced by less than 4% 
from an initial value t ps 150[ieV. 

In addition to differences in the peak shape the tem- 
perature dependence of dl/dV was also distinct in the 
two regimes (see Fig. 4). In the single peak regime, the 
zero-bias dl/dV, g(T), decreased logarithmically with T 
(Fig. 4(a)). In contrast, when double peaks appeared, 
g(T) exhibited a non-monotonic behavior (Fig. 4(b)), 
where with increasing T g(T) slightly increased initially, 
then slowly decreased before increasing again when T 
exceeded Tk- Together, these evidences point to quali- 
tatively different phases in the two regimes and the exis- 
tence of a quantum transition between them. 

To date no theoretical work has addressed the parallel- 
coupled DQD with inter-dot tunnel coupling. Never- 
theless, because the ZBA occurs under slightly non- 
equilibrium conditions it is likely that we may identify 
the observed transition with the quantum critical phe- 
nomenon discussed in the two-impurity Kondo problem 



for the series geometry based on thc following evidence: 
(a) a continuous evolution from the single- to double 
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FIG. 4: Temperature dependence of the differential conduc- 
tance dl/dV versus Vsd in the Kondo valley of Fig. 3(c): (a) 
First quantum state regime: from top to bottom T=47, 70, 
95, 110, 150, 180, 220, 280 mK, at V5=-0.5953V. (b) Second 
quantum state regime: from bottom to top T= 50, 80, 100, 
160, 180,210, 250, 270 mK, at V5=-0.5928V. Insets show the 
distinctly different temperature dependences of the zero-bias 
linear conductance, g(T) = dI/dV\y=o, in the two regimes. 



peaked behavior, (b) a maximum in g, the zero-bias 
dl/dV, near the transition point, (c) different behaviors 
in the temperature dependence of g, and (d) a strong 
renormalization of the peak splitting, 6, compared to the 
estimated t or AF coupling, J, close to the transition. 
These features are all in qualitative agreement with pre- 
dictions for a series-coupled DQD [llllllllllllllll, 
although ideally theory predicts a maximum g reaching 
the unitarity limit 4e 2 /h. Semi-quantitively, it is infor- 
mative to roughly estimate key parameters and compared 
these to observed splitting, 6 < 120/ieV. Within the the- 
oretical scenarios, S must be compared to At ss 600/ieV 
or 2 J = 8t 2 /U w ISOfieV. (Note that in the open dot 
regime U is expected to be reduced from its closed dot 
value by roughly 1/3 yielding U « lmeV.) The re- 
duction of S compared to At and 2J are in agreement with 
theory and lends further credibility to our identification 
of thc quantum transition. 
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